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ABSTRACT 

Aims. Our goal was to improve the fundamental parameters of the massive close double-lined eclipsing B2.5V+B2.5V binary 
CV Velorum. 

Methods. We gathered new high-resolution echelle spectroscopy on 13 almost consecutive nights covering essentially two orbits. 
We computed a simultaneous solution to all the available high-quality radial-velocity and light data with the latest version of the 
Wilson-Deviney code. 

Results. We obtained the following values for the physical parameters: Mi = 6.066(74)M G , M 2 = 5.972(70)M o , Ri = 4.126(24)fi , 
R2 = 3.908(27)i? o , logLj = 3.20(5)L o , and logL 2 = 3.14(5)L Q . The quoted errors contain a realistic estimate of systematic uncer- 
tainties mainly stemming from the effective temperature estimation. We derived abundances for both components and found them to 
be compatible with those of B stars in the solar neighbourhood. We discovered low-amplitude periodic line-profile variations with the 
orbital frequency for both components. Their interpretation requires new data with a longer time span. The primary rotates subsyn- 
chronously while the secondary's v sin i and radius are compatible with synchronous rotation. Finally, we provide an update of the 
empirical mass-luminosity relation for main-sequence B stars which can be used for statistical predictions of masses or luminosities. 
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1. Introduction 

Calibrating stellar evolution models requires to know the phys- 
ical parameters of stars with a very high accuracy. Eclipsing 
and spectroscopic binaries with detached components ensure to 
obtain the required precise parameters. In this respect, CV Vel 
(HD 77464, Vmag = 6.7, spectral type B2.5V+B2.5V) and a 
few other massive binary stars have an important place in astro- 
physics. The star CV Vel is one of the well-known double-lined 
early-type eclipsing binary systems. 

The binarity of CV Vel was originally established by van 
Houten (1950). The light curve of the system was first obtained 
by Gaposchkin (1955) photographically, leading to the first es- 
timate of the physical parameters of the components. The sys- 
tem was studied spectroscopically by Feast (1954) and Andersen 
(1975). Feast (1954) obtained 34 spectra and derived the first 
spectroscopic orbital parameters of the system. He reported the 
spectral type of the components as B2V+B2V. Andersen (1975) 
performed a spectral study of the system based on data spread 
over 83 days and obtained much more accurate orbital parame- 
ters. Clausen & Gr0nbech (1977) subsequently obtained a high 
quality light curve of CV Vel in the Stromgren bandpasses. 
These authors derived the orbital and physical parameters of the 
system adding the radial velocities of Andersen (1975) to their 
data. The orbital period of the star amounts to 6.889494(8) days. 

CV Vel is an important massive binary whose parameters can 
in principle be improved from modern high-quality data, which 
is why we revisited this star after 30 years. We present a new de- 
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tailed analysis of CV Vel considering all the available archival 
data to which we add new high-quality echelle spectra covering 
two consecutive weeks. We also search for a signature of oscilla- 
tions by carefully analysing the line profiles of the components, 
given that they reside in the instability strip of the slowly pulsat- 
ing B stars. 

2. New spectroscopic observations 

New spectra of CV Vel were obtained using the Swiss 1.2-m 
Euler telescope at La Silla, Chile and the CORALIE spectro- 
graph. The system was observed on 13 almost successive nights 
in December 2001 and January 2002. TableQ] lists the spectral 
lines used in our radial-velocity analysis. The resolution of the 
spectra is 50000 and the integration time was in each case one 
hour to achieve an average S/N ratio near 200 in the blue part of 
the spectrum. 

The CORALIE spectra were reduced using the standard re- 
duction pipeline (Baranne et al., 1996). For the study of the or- 
bital motion, Gaussian fits to the absorption lines were made in 
order to derive their central wavelengths, thus transforming into 
the radial velocities for each of the considered lines. To illustrate 
the data quality, we show in Fig.[T]a selection of four lines for 
the two orbital phases with the largest line separation. It is clear 
from this figure that the primary has narrower lines than the sec- 
ondary while being of similar temperature. This means that one 
of them must be asynchronous if the radii are similar, as we 
show is the case below. We come back to this point in Sect. 5. 

We list in Table|2] the average value of the measured radial 
velocities obtained from the five least blended spectral lines, i.e. 
He I 4120.993A, Mg II 4481.327 A, Si III 4552.622A, Si III 
4567.840A and He I 6678. 149 A lines. For each of these lines, 
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Table 2. Radial velocity values for CV Vel computed in this 
study. 
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Fig. 1. Observed Mgll 4481, Sim 4552, Ha and Hel 6678 line 
profiles at orbital phases 0.29 and 0.75. 

Table 1. Wavelength ranges and the spectral lines considered in 
our study. 
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the orbit was solved separately and the obtained results are sum- 
marized in Table[3] The internal precision of a single observa- 
tion amounts to ~ 1.0 kms" 1 for the primary and ~ 2.0kms~ 1 
for the secondary. We note a quite significant deviation among 
the Vo values derived from the different lines. This is not un- 
common in this part of the HR diagram (e.g. Mathias et al. 
(1994) for a discussion). It is interpreted in various ways, among 
which unknown blends within the lines, an asymmetric Stark ef- 
fect among different lines (Struve & Zebergs 1960) or the pres- 
ence of a so-called Van Hoof effect (Van Hoof 1975) in pulsat- 
ing stars. The latter effect is due to the fact that spectral lines 
form in different line forming regions in the stellar atmosphere, 
such that a running wave nature of the oscillations implies a shift 
in phase during the velocity cycle. Such an effect has clearly 
been detected in the pulsating B stars B W Vulpeculae and a Lupi 
(Mathias & Gillet 1993). Our data have insufficient time base 
and temporal resolution to test the presence of a Van Hoof ef- 
fect in CV Vel (see also Section 6). Even if a Van Hoof effect is 
present, it would not explain the difference in Vo for the two Si 
lines of the same triplet (Silll /14553 and Silll ,14568). A similar 
unexplained average velocity shift between the same lines in this 
Si triplet was encountered by De Cat et al. (2000) in their analy- 
sis of several hundred line-profile data, of very high S/N, of the 
B0.5III star s Per. As there, we are unable to offer a good inter- 
pretation of this shift. We will take it into account as a systematic 
uncertainty in the derivation of the orbital and fundamental pa- 
rameters of the star (see below). Though the eccentricity, e, was 
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first regarded as a free parameter, its value always appeared to 
be zero, within the error range, so we do not list it in Table[3] 

The fifth column of Table|4] gives the orbital solution from 
the averaged radial velocities of our study. The average radial- 
velocity values are plotted versus the orbital phase for each of 
the components in Fig. [2] The phases were calculated using the 
ephemeris given by Clausen & Gr0nbech (1977): 

Pri. Min. = HJD 2442048 d 66894(14) + 6^889494(8) x E. (1) 



3. Simultaneous radial-velocity and light curve 
analysis 

To measure the £ a>(0 - C) 2 values, given in Table|4j we have 
re-analyzed the data of Feast (1954) and Andersen (1975). This 
showed that the data of Andersen (1975) and those of our pa- 
per have similar quality while those of Feast (1954) have clearly 
lower quality, as already emphasised by Andersen (1975). For 
this reason, we omitted Feast's data in the remaining of this pa- 
per. 

We combined the newly obtained radial-velocity values 
(Table 2) with those previously obtained by Anderson (1975). 
This gave a total of 62 radial-velocity values for each of the com- 
ponents to which we assigned equal weights. As can be seen in 
TableHJ we find an offset in average velocity Vo between the two 
radial- velocity sets. This could be due to the difference in used 
spectral lines in the two studies, but it may also be that we find 
a downward trend in Vo over time due to the presence of a third 
body. We have too few epochs to model any Vq trend for the time 
being. For this reason, we shifted the data towards the CORALIE 



K. Yakut, C. Aerts and T. Morel: The early-type close binary CV Velorum revisited 3 

Table 3. Spectroscopic orbital elements of CV Vel for different spectral lines. See text for details. The formal standard errors cr are 
given in parentheses in the last digit quoted. 
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Table 4. Spectroscopic orbital parameters of CV Vel. The standard errors cr are given in parentheses in the last digit quoted. 
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Fig. 2. The radial velocity observations of CV Vel versus the or- 
bital phase. The filled and open circles correspond to the veloci- 
ties of the primary and the secondary, respectively. The data are 
listed in Tableland the sine curve corresponds to the elements 
given in Table|4] The center-of-mass velocity is indicated by the 
dashed-dotted line. 



Vq value, keeping in mind that this systematic uncertainty adds 
to the uncertainty of the overall solution. 

These 62 radial-velocity values were then combined with the 
light curves obtained in the Stromgren (ubvy) filters by Clausen 
& Gr0nbech (1977), kindly made available to us by these au- 
thors. The u band data show a somewhat larger scatter than the 
other filters so we used different weights. We determined the 
weights for the uvby bands from the observing errors given by 
Clausen & Gr0nbech (1977), resulting in the values 1.7, 2.5, 2.0, 
and 2.0 for u, v, b and y, respectively. 

The previous light curve analyses of the system have been 
done using the WINK method by Clausen & Gr0nbech (1977) 



(second column of Table 5). Clausen & Gr0nbech (1977) have 
used a linear limb-darkening law, fixing the temperature of the 
primary star at the value 18300 K. They determined solutions 
for each colour, assuming different values for k — ^ (see, their 
Table 4). In this work, the light and radial-velocity curves were 
analyzed simultaneously using the 2005 September version of 
the Wilson-Deviney code (hereafter WD, Wilson & Devinney 
1971; Wilson 1994). Mode 2 of the WD code, which assumes 
that both of the components are detached, was adopted. In the 
final solution, we fixed the following parameters: the limb- 
darkening coefficients (taken from Diaz-Cordoves & Gimenez 
1992, see Table 5), the values of the gravity darkening coeffi- 
cients and the albedos. The adjustable parameters are the incli- 
nation i, the temperature of the secondary component T%, the 
luminosities L\b and L\ y , and the surface potentials Q.\ and Q2. 
The differential correction (DC) code was ran until the correc- 
tions to the input parameters were lower than their errors. 

To obtain a simultaneous solution, different possibilities 
were tested. We varied the primary's temperature in the range of 
17800 K - 19000 K with 50 K steps and found that Ti = 18000 K 
gave the best solution with the lowest formal errors for the un- 
known quantities and the lowest £(0 - C) 2 . However, the solu- 
tions with T\ between 18000 K and 19000 K were almost indist- 
inghuishable. This is in agreement with the solution proposed by 
Clausen & Gr0nbech (1977), who assigned T x = 18200 K to the 
mean component of the system and considered an uncertainty of 
500 K. In selecting the best solution, we first assumed the sys- 
tem to be eccentric leaving the e and a> values free. As a result 
we obtained always e < 0.001. The final results are given in 
Table|5]for each band separately, as well as for the simultaneous 
solution including all photometric data. We find a similar orbital 
inclination value than Clausen & Gr0nbech (1977). The com- 
puted light and radial-velocity curves for the parameters given 
in Table|5] are shown by the solid lines, and are compared with 
all the observations, in Fig. [3] 

4. Abundance determination 

Our spectra are of sufficient quality, and cover a wide enough 
wavelength range, to undertake an abundance analysis for the 
two components of the CV Vel system. The non-local thermo- 



4 



K. Yakut, C. Aerts and T. Morel: The early-type close binary CV Velorum revisited 



Table 5. Orbital parameters and photometric elements, with their formal errors, of CV Vel. See text for details. CG: Clausen & 
Gr0nbech (1977), the other columns represent results obtained in this work. The parameters with a superscript * were kept fixed. 
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0.9835(10) 


Fractional radii of the primary 














•\ pole 


- 


0.1172(4) 


0.1180(3) 


0.1181(3) 


0.1182(3) 


0.1172(1) 


' 1 point 


- 


0.1177(4) 


0.1186(3) 


0.1186(3) 


0.1187(3) 


0.1184(1) 


M side 


- 


0.1174(4) 
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0.1122(3) 
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0.113(1) 
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Radiative parameters: 
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0.028 


0.026 
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dynamic equilibrium (NLTE) abundances of He, C, N, O, Mg, 
Al, Si, S and Fe were calculated using the latest versions of 
the line formation codes DETAIL/SURFACE and plane-parallel, 
fully line-blanketed Kurucz atmospheric models (Kurucz fl993l l. 
Curve-of-growth techniques were used to determine the abun- 
dances using the equivalent widths (EWs) of a set of spectral 
lines measured on a mean spectrum created by averaging the five 
CORALIE spectra taken at 0-0.73 (see Morel et al. l2006l for the 
line list used). These EWs were subsequently corrected to ac- 
count for the dilution by the companion's continuum, assuming 
the luminosity ratio given in Table|7] The two components have 



essentially the same effective temperature, so that these correct- 
ing factors can be considered independent of wavelength. 

We first computed the abundances assuming the logg and 
r e ff values obtained from the light curve analysis (Table|7J, but 
found that the Si ii/m ionization balance is not well fullfilled in 
that case. Instead, a slightly higher value of T e g = 19 000 ± 500 K 
is needed for both components. Because of this ambiguity sur- 
rounding the temperature determination, we have estimated the 
abundances using both the photometric and the spectroscopic 
values (hereafter r p h to and ^ionization, respectively). The uncer- 
tainties on the elemental abundances were calculated by adding 
in quadrature the internal errors (i.e. the line-to-line scatter) and 
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u + 0.25 




(c) 

Fig. 3. The observed and computed (a) light and (b) radial- 
velocity curves of CV Vel. The curve in the u, b and y bands are 
moved by +0.25, -0.25 and -0.50, respectively, in intensity for 
good visibility. Circles and squares denote the data in this study, 
and Anderson (1975) respectively. The center-of-mass velocity 
is indicated by the dashed-dotted line, (c) 3-D representation of 
CV Vel at phases 0.00, 0.25 and 0.50. 



the errors arising from the uncertainty on the microturbulent ve- 
locity, The latter quantity is usually inferred in B stars by 
requiring no dependence between the abundances yielded by 
lines of an ionic species with strong features (e.g. On) and their 
strength. However, most metallic lines are intrinsically weak in 
the relevant temperature range (measured EWs < 20 mA), such 
that it is virtually impossible to constrain the microturbulence. 
We adopt in the following a canonical value £=5+5 km s , with 
the large conservative error bars encompassing all plausible val- 
ues for B dwarfs. This translates into poorly defined abundances 
for the elements solely exhibiting relatively strong lines, in par- 
ticular He, Si and Mg. 

The results for the two choices of r e ff are provided in Table[6] 
No abundances are found to deviate at more than the 3cr level 
from the typical values found for early B dwarfs in the solar 
neighbourhood (Daflon & Cunha 2004) for both T e g estimates. 



Therefore, there is no convincing evidence for departures from 
a scaled solar composition. The nominal abundance values are 
systematically higher in the secondary than in the primary. As 
the chemical composition of the two components is expected 
to be identical for such a young, detached system with no past 
or ongoing episodes of mass transfer, this strongly argues for a 
higher microturbulent velocity in the secondary. This difference 
is, however, not sufficient to explain the discrepant rotation rates 
of the two components (see Section 6 below). 

Recent works point to the existence of a population of 
slowly-rotating B-type (magnetic) pulsators exhibiting an un- 
expected nitrogen excess at their surface, which could possi- 
bly be attributed to deep mixing or diffusion effects (Briquet & 
Morel 2007; Morel et al. 120061 1. A similar phenomenon cannot 
be firmly established in either of the two components of CV Vel: 
while the ratio of the abundances of N and C ([N/C]) is compa- 
rable to the values previously reported for the N-rich /3 Cephei 
and slowly pulsating B stars, the [N/O] ratio is indistinguishable 
from solar. This conclusion holds irrespective of the temperature 
scale adopted. 

Regarding the derivation of the physical parameters of the 
component stars, we are thus left with the conclusion that a sec- 
ond systematic uncertainty occurs, besides the one induced by 
the offset in Vq. We forced the effective temperature of the pri- 
mary to 19000 K and checked the orbital solutions from WD 
discussed in the previous section. As already mentioned, the dif- 
ferences in the orbital parameter values are small and compatible 
with the uncertainties listed there. While deriving the physical 
parameters of the system in the next section, we propagated the 
errors caused by the systematic uncertainty in effective temper- 
ature in order to achieve realistic errors. 



5. Physical parameters of the system 

The physical parameters of C V Vel first obtained by Gaposchkin 
(1955) were Mi = 5.618 M Q , M 2 = 5.418 M Q , Ri = 4.26 R Q , and 
R2 = 4.16^0. Andersen (1975) reported the parameters of the 
components as M\ 2 = 6.05 M , R\ 2 = 4.05 R Q , and log T e ff = 
4.26. Clausen & Gr0nbech (1977) subsequently obtained Mi = 
6.10 M , M 2 = 5.99 M @ , Ri = 4.10 R Q , R 2 = 3.95 fl Q , T x = 
1 8200 K and T 2 = 18060K. 

We summarize the obtained parameters of the system and of 
the components obtained from our simultaneous fit in Table|7] 
In computing them, we took the values of the light curve esti- 
mates for the effective temperature of the components. We argue 
that the photometric data are of high quality while the spectro- 
scopic temperature estimates suffer from the unknown microtur- 
bulence. However, we take into account a systematic uncertainty 
of 500 K for the temperature. This equals the uncertainty derived 
from the silicon ionisation balance and is a factor ten higher than 
the formal error that resulted from our WD analysis, so we con- 
sider this to be a safe conservative approach. We propagated this 
uncertainty for all the derived quantities relying on the effective 
temperature. 

We derived the component's bolometric corrections from the 
temperatures using the interpolation formula given by Balona 
(1994) for G5 to early-type stars. In order to find the distance of 
the system, we used the E(B-V) value of 0^08 from Savage et al. 
(1985) resulting in d — 553 pc. 

Our results are in very good agreement with those obtained 
by Clausen & Gr0nbech (1977). We were unable to achieve more 
precise values than they did, due to the uncertainty on the effec- 
tive temperature resulting from our spectroscopic analysis. 
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Table 6. Mean NLTE abundances of the primary and secondary components of CV Vel (on the scale in which log e[H] = 12), along 
with the total l-cr uncertainties. The number of used lines is given in brackets. We quote the values obtained when using the 
temperatures derived from the light curve analysis (T p hoto) or from the ionization balance of silicon (Ionization). The gravity is fixed 
in both cases to the values quoted in Tableland the microturbulent velocity, to 5 km s . For comparison purposes, the first 
column gives the typical values found for nearby B dwarfs (Daflon & Cunha 2004). We define [N/C] and [N/O] as log[e(N)/e(C)] 
and log[e(N)/e(0)], respectively. 







left — Iphoto ' 


= 18000 K 


^eff — ^ionization 


=19000 K 




OB stars 


Primary 


Secondary 


Primary 


Secondary 


He/H 


0.085 


0.053+0.029 (6) 


0.097+0.033 (6) 


0.040±0.022 (6) 


0.073+0.031 (6) 


loge(C) 


-8.2 


8.08+0.10(4) 


8.29±0.17 (2) 


7.89±0.07 (4) 


8.11±0.17 (2) 


loge(N) 


-7.6 


7.93±0.13 (7) 


8.16±0.16(6) 


7.69±0.13 (7) 


7.92±0.17(6) 


loge(O) 


-8.5 


8.74±0.18 (8) 


8.85+0.20 (4) 


8.44±0.17 (8) 


8.55±0.19(4) 


l0g6(Mg) 


-7.4 


7.01±0.40 (1) 


7.26+0.45 (1) 


7.13+0.40(1) 


7.39±0.39 (1) 


loge(Al) 


-6.1 


6.32±0.18 (3) 


6.48+0.17(3) 


6.12+0.18 (3) 


6.27±0.17(3) 


loge(Si) 


-7.2 


6.97±0.44 (7) 


7.15±0.61 (7) 


6.93±0.28 (7) 


7.11+0.44(7) 


loge(S) 


-7.2 


7.02±0.21 (5) 


7.25±0.23 (2) 


7.08±0.19(5) 


7.29±0.25 (2) 


log e(Fe) 


-7.3" 


7.20+0.20 (3) 


7.46+0.18(2) 


7.05±0.15 (3) 


7.34±0.23 (2) 


[N/C] 


—0.6 


-0.15+0.17 


-0.13±0.24 


-0.20±0.15 


-0.19±0.25 


[N/O] 


—0.9 


-0.81±0.23 


-0.69±0.26 


-0.75±0.22 


-O.63±0.26 



° From Morel et al. ( 2006 1. 



Table 7. Absolute parameters for CV Vel. The errors are given in parentheses and take into account the systematic uncertainties 
encountered for V and for the effective temperature of the components. The effective temperature of the Sun was taken as 5780 K, 
its bolometric absolute magnitude as 4.75 mag and its bolometric correction as -0.07 mag. 





Unit 


primary component 


secondary component 


Mass (M) 


M 


6.066 (74) 


5.972 (70) 


Radius (R) 


Ro 


4.126 (24) 


3.908 (27) 


Effective temperature (log r cff ) 


K 


4.255(12) 


4.250(12) 


Luminosity (log L) 


I- 


3.204(48) 


3.137(48) 


Surface gravity (log g) 


cgs 


3.99 (6) 


4.03 (6) 


Bolometric correction (BC) 


mag 


-1.68 (5) 


-1.66 (5) 


Absolute bolometric magnitude (M&oi) 


mag 


-3.26(3) 


-3.09 (3) 


Absolute visual magnitude (My) 


mag 


-1.58(6) 


-1.43 (5) 


Semi-major axis (a) 


Ro 


34.90(15) 




Distance (d) 


pc 


553 (32) 





We confronted our observational results for CV Vel with evo- 
lutionary models computed by De Cat et al. (2006) with the Code 
Liegeois d'Evolution Stellaire (CLES, Scufiaire 2005). These 
computations were done using the new solar abundances as re- 
ported by Asplund et al. (2005) and a standard mixing length 
description of convection with l m = 1 .75 times the local pres- 
sure scale height. The use of these models is appropriate for 
the components of CV Vel. Indeed, the Z-values we obtained 
for the primary and secondary are Z = 0.0119 + 0.0028 and 
Z = 0.0160 + 0.0040 for the case of an effective temperature 
of 18 000K and where we have taken the standard abundances 
of Grevesse & Sauval (1998) for the elements not included in 
Table 6. The metallicities we get are thus consistent with those 
used in the model computations by De Cat et al. (2006). The 
tracks for M = 5.9, 6.0, 6.1 M are shown in Fig.|4]for two cases: 
no core overshoot and an overshoot of 0.2 times the local pres- 
sure scale height. We overplotted the position of CV Vel's com- 
ponents according to our results mentioned in Table[7] It can be 
seen that the agreement between our observational results and 
the evolutionary models is satisfactory. Recent results from as- 
teroseismology have shown that a small core overshoot value is 
necessary to bring observed and identified oscillation frequen- 
cies in agreement with evolutionary models (Aerts et al. 2003; 
Pamyatnykh et al. 2004) in this part of the HR diagram. The sys- 
tematic uncertainty for the temperatures (and thus luminosities) 



prevented us from constraining the overshoot parameters for the 
components of CV Vel. We derive an age of about 40 million 
years for the system. Clausen & Gr0nbech (1977) reported an 
age of only about 30 million years, but they used of course much 
older models with somewhat different input physics (e.g., older 
opacity values). 



6. Study of the line-profile variability 

The components of CV Vel reside in the instability strip of the 
slowly pulsating B stars, not far away from the /JCephei strip 
(Fig. [4]). It is therefore worthwhile to search for short-period 
spectroscopic variations of the components. Percy & Au- Young 
(2000) reported the possibility of a variation with an amplitude 
of 0.02 mag and a period of 0.25 day in the Hipparcos photom- 
etry. However, they did not study CV Vel in detail. Clausen (pri- 
vate communication) was unable to find short-period variations 
in his extensive photometric data. 

In Fig.|5^ we plot the Mg II 4481 A lines of the components 
obtained at the phases outside of the eclipses. This line is the 
best one to study line-profile variability because it is one of the 
deeper ones and has the best S/N level. We shifted all the lines to 
the common center of mass. The bold continuous line in Fig.[5t> 
is the average profile. In Fig. [5]:, the residuals of each line with 
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Fig. 4. The components of CV Vel in the HR diagram with error 
bars. The upper point represents the primary with M = 6.07 M 
and the lower point indicates the secondary with M = 5.97 M 
(see TableO. The full thick line indicates the /3 Cephei instability 
strip, the full dashed line is the instability strip of the slowly pul- 
sating B stars. Both strips were computed with the non-adiabatic 
oscillation code MAD (Dupret, 2001) assuming that no core 
overshoot occurs. Evolutionary tracks computed with CLES are 
shown for M = 5.9, 6.0, 6.1 M & , and for zero core overshoot 
(thin full lines) as well as a core overshoot of 0.2 (thin dashed 
lines). 




1.00 

0.90 



4481.0 4481.5 



3.5 4481.0 4481 



482.0 4482 




Fig. 5. (a) Normalized Mg II(/i4481) line profiles of the primary 
(left panel) and the secondary (right panel). Phase runs from the 
top to the bottom, (b) Line profiles superimposed. The full thick 
line is the average profile, (c) Residual profiles with respect to 
the average 



respect to the average are shown. It is apparent from this figure 
that the primary shows line-profile variations. 

We first of all estimated the rotational and thermal broad- 
ening from the average profiles shown in Fig.[5]3. It is evident 
from this figure, as well as from Fig.Q] that the secondary rotates 
faster than the primary since both components essentially have 
the same thermal broadening. This is in contrast to Andersen's 
(1975) conclusion that the components have equal rotation ve- 
locity. We computed theoretical Mg II profiles for Gaussian in- 
trinsic broadening as well as rotational broadening and find 
vsinz = 19 + lkms" 1 for the primary and 31+2kms _I for 
the secondary. Assuming the orbital and rotational axes to be 
aligned, and using the radius estimates from Table 6, we find a 
rotation period of 11.0 (+0.6) days for the primary and of 6.4 
(+0.4) days for the secondary. It is noteworthy that the primary 
rotates subsynchronously, while the secondary's rotation period 
is compatible with the orbital period within the errors. Tidal 
theory (e.g., Zahn 1975, 1977) predicts the occurrence of syn- 
chronisation before circularisation, and both to happen on time 
scales shorter than about half the main-sequence duration, de- 
pending on the birth values of the eccentricity and orbital pe- 
riod. Deviations from synchronisation in circularised binaries 
have been reported before for close binaries with a B-type com- 



ponent and a similar orbital period, e.g. in the binary ^Ori with 
a pulsating component (De Mey et al. 1996). 

We tried to quantify the line-profile variability of both com- 
ponents by computing the line moment variations in the defini- 
tion of Aerts et al. (1992) for the 23 spectra shown in Fig.|5h. 
These line diagnostics stand for the centroid variation (< v >), 
the variation in the line width (< v 2 >) and the variation in the 
line skewness (< v 3 >). They are specifically powerful to unravel 
low-degree (I < 4) non-radial oscillation modes. We performed 
a frequency search with the Scargle method (Scargle 1982) on 

< v > in the range [0.0, 3.4] c d _1 , where the upper limit of this 
interval corresponds to the Nyquist frequency. As a result, our 
dataset is not suitable to detect p-mode oscillations in B2.5V 
stars because those have frequencies of typically 4 to 8 c cT 1 for 
slow rotators (e.g. Aerts & De Cat 2003). We should be able to 
find a signature of g-modes with amplitudes of a few kms~' as 
in slowly pulsating B stars (e.g. De Cat & Aerts 2002). 

The highest frequency peak in the Scargle periodogram of 

< v > for the primary occurs at / = 0. 166+0. 009c d~', with an 
amplitude of 2.2+0.4 km s _I . This frequency has a variance re- 
duction of 54%, which implies a decrease in the standard devia- 
tion from 2.3 to 1.5kms _I . It reaches a level of 3.8 S/N, where 
the S/N level was computed from an oversampled Scargle peri- 
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Table 8. Observed parameters for B-type main-sequence double-lined eclipsing binary stars. 



Star 


Sp.T 


P 


Mi 


M 2 


Ri 


^2 


log T e \ 


log T e 2 


logL! 


log Li 


Ref 


AH Cep 


B0.5Vn+B0.5Vn 


1.78 


15.4 


13.6 


6.38 


5.86 


4.476 


4.456 


4.465 


4.314 


1 


CW Cep 


B0.5V+B0.5V 


2.73 


13.484 


12.05 


5.685 


5.177 


4.452 


4.442 


4.27 


4.15 


2 


QXCar 


B2V+B2V 


4.48 


9.245 


8.46 


4.289 


4.051 


4.377 


4.354 


3.72 


3.58 


2 


V497 Cep 


B3V+B3V 


1.20 


6.89 


5.39 


3.69 


2.92 


4.290 


4.249 


3.265 


2.860 


3 


V539 Ara 


B3V+B4V 


3.17 


6.239 


5.313 


4.432 


3.734 


4.26 


4.23 


3.29 


3.02 


2 


CV Vel 


B2.5V+B2.5V 


6.89 


6.066 


5.972 


4.126 


3.908 


4.255 


4.250 


3.20 


3.14 


4 


AG Per 


B4V+B5V 


2.03 


5.36 


4.9 


2.99 


4.297 


4.260 


4.241 


2.95 


2.75 


5 


UOph 


B5V+B5V 


1.68 


5.186 


4.672 


3.438 


3.005 


4.22 


4.199 


2.91 


2.7 


2 


DI Her 


B4V+B5V 


10.55 


5.173 


4.523 


2.68 


2.477 


4.23 


4.179 


2.73 


2.46 


2 


V760 Sco 


B4V+B4V 


1.73 


4.968 


4.609 


3.013 


2.64 


4.228 


4.21 


2.82 


2.63 


2 


GGLup 


B7V+B9V 


1.85 


4.155 


2.532 


2.411 


1.753 


4.176 


4.041 


2.42 


1.61 


2 


( Phe 

* 2 Hya 


B6V+B8V 


1.67 


3.92 


2.545 


2.851 


1.853 


4.163 


4.076 


2.51 


1.79 


2 


B8V+B8V 


2.27 


3.605 


2.631 


4.384 


2.165 


4.066 


4.041 


2.5 


1.79 


2 


IQ Per 


B8V+B6V 


1.74 


3.513 


1.732 


2.446 


1.503 


4.09 


3.885 


2.09 


0.85 


2 


PVCas 


B9.5V+B9.5V 


1.75 


2.82 


2.761 


2.243 


2.287 


4.000 


4.000 


1.65 


1.67 


2 


V451 Oph 


B9V+A0V 


2.20 


2.769 


2.36 


2.64 


2.03 


4.033 


3.991 


1.93 


1.53 


6 


GG Ori 


B9.5V+B9.5V 


6.63 


2.342 


2.338 


1.852 


1.830 


3.9978 


3.9978 


1.480 


1.470 
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References for Table 7.: 1 Holmgren et al. (1990); 2 Andersen (1991) and references therein ; 3 Yakut et al. (2003); 4 This study; 5 Gimenez & 
Clausen (1994); 6 Clausen et al. (1986); 7 Torres et al. (2000) 
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Fig. 6. The first moment variation as a function of orbital phase 
for the primary (o) and secondary (A). The full lines with corre- 
sponding symbols are a least-squares sine fit. 

odogram of the residuals over the entire interval [0.0, 3.4] c d -1 . 
We recover the same frequency / in < v 3 >. We do not, how- 
ever, find any periodic phenomenon in < v 2 >. A pulsation mode 
would imply the occurrence of / and/or 2f in < v 2 > as well 
(Aerts et al. 1992, De Cat et al. 2005). It may be that our dataset 
is too limited in quality to detect a frequency in < v 2 > given that 
the even moments are always noisier than the odd ones (Aerts et 
al. 1992). 

For the secondary, we find /' = 1.168+0. 009 cd~'. This 
is a one-day alias of /. Forcing / leads to an amplitude of 
1.8±0.4kms _I (at level 4.1 S/N) and a variance reduction of 
51%, bringing the standard deviation from 1.5 to 0.9kms _I . 
This time, we cannot recover any frequencies in < v 2 > or 
< v 3 >. We show in Fig. [6] the centroid velocity < v > for both 
the primary and secondary, folded with the frequency /, as well 
as their least-squares fit. We find that the primary's < v > lags 
behind a quarter of a phase compared to the one of the secondary. 

We conclude to have found evidence of periodic line-profile 
variability in the primary and secondary of CV Vel with a period 



between 5.7 and 6.4 days. This is almost equal to the orbital pe- 
riod. Moreover, the quoted frequency error is a formal lcr least- 
squares error, and is an underestimation of the true error given 
that we covered less than two full cycles with unblended Mg II 
lines. It therefore seems that line-profile variability occurs with 
a period very close or equal to the orbital one. This is potentially 
interesting as it could point towards the excitation of a tidally 
induced g-mode oscillation in the components. However, given 
the phase relation we obtained, and the absence of any periodic 
signature in < v 2 >, it is more likely that a reflection and/or 
variable limb-darkening effect lies at the origin of the detected 
variability. We need a more extensive data set to firmly establish 
the correct interpretation of CV Vel's line-profile variability. 

7. Summary 

From combined existing light and old and new radial-velocity 
curves, we have obtained a full solution of the double-lined 
eclipsing binary CV Vel. The system's center-of-mass velocity 
(Vo) obtained by Feast (1954), Andersen (1975) and in this study 
are 27.9(1.1), 24.3(5), and 21.7(2) km s -1 , respectively. These 
three different and decreasing values may indicate the presence 
of a distant third body orbiting the binary system. Forthcoming 
observations of the system are necessary to evaluate this hypoth- 
esis. We checked but did not find a significant orbital period 
change. 

We have corrected for the offset in Vq in Andersen's (1975) 
and our data and merged these two high-quality radial-velocity 
curves to improve our knowledge on the fundamental param- 
eters of the components. In order to achieve this, we comple- 
mented the spectroscopic data with the high-quality 4-colour 
(uvby) lightcurves obtained by Clausen & Gr0nbech (1977). 
These photometric and radial-velocity data were then solved si- 
multaneously. The eclipsing and detached properties of the sys- 
tem have led to reliable orbital and physical parameters in agree- 
ment with the values obtained earlier by Clausen & Gr0nbech 
(1977). 

We presented for the first time an abundance analysis for the 
system and found results in agreement with those for B stars 
in the solar neighbourhood for both components. This analy- 
sis led to a systematic uncertainty of some 500 K for the effec- 



K. Yakut, C. Aerts and T. Morel: The early-type close binary CV Velorum revisited 



9 



tive temperatures of the components. Thus, these two parame- 
ters of this seemingly well-known system remain uncertain de- 
spite our detailed analysis based on echelle spectroscopy. It is 
not uncommon to find systematic uncertainties and deviations 
between photometrically and spectroscopically derived effective 
temperatures of B-type stars (e.g. De Ridder et al. 2004 for a 
discussion), and even to have uncertainties of order 500 K for 
such stars among methods based on the same data (e.g. Smalley 
& Dworetsky 1995, Morel et al. 2006, Kaiser 2006). Previous 
works on CV Vel did not include a spectroscopic temperature 
estimate and were thus not able to estimate the systematic un- 
certainty of the component's temperatures. 

Finally, we discovered line-profile variability in both com- 
ponents. The variability is most prominent in the primary. From 
the characteristics of the line profile moments, we tentatively 
interpreted this variability as due to an extrinsic orbital effect 
rather than an intrinsic oscillation. Any final interpretation of the 
detected line-profile variability requires a much more extensive 
dataset, however. 

Our results show that CV Vel consists of young stars which 
are about half-way in their central hydrogen burning phase. 
Though the orbit is circular, the components clearly have dif- 
ferent rotation speeds. We thus conclude that the tidal forces had 
insufficient time so far to bring the system to a circular orbit with 
co-rotating components, i.e. the circularisation process is com- 
pleted but not the synchronisation one. 

Observational tests of stellar evolution models have been 
done by many authors previously (Pols et al. 1997, Young 
et al. 2001, Lastennet & Valls-Gabaus 2002, Young & Arnett 
2005, and references therein). To determine the masses and 
radii of stars independently from models, one should use de- 
tached, eclipsing double-lined binaries and compute a simul- 
taneous solution to their radial-velocity and multi-colour curves. 
Unfortunately, accurate physical parameters are only available 
for very few systems with early-type components. Andersen 
(1991) compiled a list of stars whose physical parameters are 
well defined, with relative mass and radius uncertainties less 
than 2%. We collected all the well-known detached eclipsing bi- 
naries with main-sequence B-type components since then from 
the literature (Tabled and added our new estimates of CV Vel. 
While we obtain a mass estimate with a relative precision of 
1.2%, the systematic uncertainty for the effective temperature 
that we derived from our high-quality spectra implies an un- 
certainty of some 12% on the luminosities of the components. 
Nevertheless, the values we obtained for CV Vel are fully com- 
patible with evolutionary models and with those of similar sys- 
tems whose parameters are known with better precision. 

The stars given in Table[8] are drawn in a mass-luminosity 
plot in Fig. |7] From Table[8]we find the following update of the 
mass-luminosity relation for the main-sequence B-type stars: 

log(L/L ) = 3.724(80) log(M/M ) + 0.162(57). (2) 

Previously, Hilditch & Bell (1987) presented a similar relation 
for O and B systems. The rms error in log(L) they found was 
0.17 whereas it is only 0.11 for our Eq. because we con- 
sider only binaries with B star components here. This equation 
is useful to predict masses from luminosities, or vice versa, for 
main-sequence B stars. 
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